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f**^ ^ The nucleon has been used as a laboratory to investigate its own spin structure and 

■""^ , Quantum Chromodynamics. New experimental data on nucleon spin structure at low to 

^ ' intermediate momentum transfers combined with existing high momentum transfer data 

offer a comprehensive picture of the transition region from the confinement regime of the 
theory to its asymptotic freedom regime. Insight for some aspects of the theory is gained 
o : by exploring lower moments of spin structure functions and their corresponding sum rules 

^ ■ (i.e. the Gerasimov-Drell-Hearn, Bjorken and Burkhardt-Cottingham). These moments 

are expressed in terms of an operator-product expansion using quark and gluon degrees 
of freedom at moderately large momentum transfers. The sum rules are verified to good 
accuracy assuming that no singular behavior of the structure functions is present at very 
1 high excitation energies. The higher-twist contributions have been examined through 

5^ ' the moments evolution as the momentum transfer varies from higher to lower values. 

' Furthermore, QCD-inspired low-energy effective theories, which explicitly include chiral 

symmetry breaking, are tested at low momentum transfers. The validity of these theories 
is further examined as the momentum transfer increases to moderate values. It is found 
that chiral perturbation calculations agree reasonably well with the first moment of the 
spin structure function gi at momentum transfer of 0.1 GeV'^ but fail to reproduce the 
neutron data in the case of the generalized polarizability 5lt- 
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1. Introduction 

In the last twenty-five years the spin structure of the nucleon led to a very produc- 
tive experimental and theoretical activity with exciting results and new challengeJ^^. 
This investigation has included a variety of aspects, such as testing Quantum Chro- 
modynamics (QCD), the theory of strong interactions, in its perturbative regime 
?;ia spin sum rules (like the Bjorken sum rule ) and understanding how the spin of 
the nucleon is built from the intrinsic degrees of freedom of the theory, quarks and 
gluons. Recently, results from a new generation of experiments performed at Jeffer- 
son Lab seeking to probe the theory in its non-perturbative and transition regimes 
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have reached a mature state. The low momentum-transfer results offer insight in a 
region known for the collective behavior of the nucleon constituents and their in- 
teractions. In this region it has been more economical to describe the nucleon using 
effective degrees of freedom like mesons and constituent quarks rather than current 
quarks and gluons. Furthermore, distinct features seen in the nucleon response to 
the electromagnetic probe, depending on the resolution of the probe, point clearly 
to different regimes of description, i.e. a scaling regime where quark-gluon corre- 
lations are suppressed versus a coherent regime where long-range interactions give 
rise to the static properties of the nucleon. ^^^^^^^^^^^ 

In this review we describe an investigatioiPlEIlZlHMMniof the spin structure 
of the nucleon through the measurement of the helicity-dependent photoabsorption 
cross sections or asymmetries using virtual photons across a wide resolution spec- 
trum. These observables are used to extract the spin structure functions gi and 52 
and evaluate their moments. These moments are powerful tools to test QCD sum 
rules and unravel some aspects of the quark-gluon structure of the nucleon. 



2. Sum rules and Moments 

Sum rules involving the spin structure of the nucleon offer an important opportunity 
to study QCD. In recent years the Bjorken sum rule at large (4-m qmentu m 
transfer squared) and the Gerasimov, Drell and Hearn (GDH) sum rul J^^ I ^'^ at 

= have attracted large experimental-'^'^ and theoretical^^ efforts that 

have provided us with rich information. This first type of sum rules relates the 
moments of the spin structure functions (or, equivalently, the spin-dependent total 
photoabsorption cross sections) to the nucleon's sta tic pro perties. The second type 
of s um rul es, such as the generalized GDH sum rul j^^ ^ ^^ or the polarizability sum 
rulePEB, relate the moments of the spin structure functions to real or virtual 
Compton amplitudes, which can be calculated theoretically. Both types of sum 
rules are based on "unsubtracted" dispersion relations and the optical theoreni^Sl. 
The first type of sum rules uses one more general assumption, such as a low-energy 
theoren|22l for the GDH sum rule and Operator Production Expansion (OPE)^ 
for the Bjorken sum rule, to relate the Compton amplitude to a static property. 

The formulation below follows closely Ref. . Consider the forward doubly- 
virtual Compton scattering (VVCS) of a virtual photon with space-like four- 
momentum q, i.e., — — (f' = < 0, where v is the energy and g is the three 
momentum of the virtual photon. The absorption of a virtual photon on a nucleon is 
related to inclusive electron scattering. The inclusive cross section, assuming parity 
conservation, contains four partial cross sections (or structure functions): ctt, ctl, 
o'TT, cr^T, (or Fi, F2, gi, 52 )• The first two are spin-averaged, while the last two 
are spin-dependent. In this review, we will concentrate on the spin-dependent ones. 
In the following discussion, we will start with the general situation, i.e., sum rules 
valid for all Q^, then discuss the two limiting cases at low and at high Q^. 

Considering the spin-fiip VVCS amplitude gxT and assuming it has an appropri- 
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ate convergence behavior at high energy, an unsubtracted dispersion relation leads 
to the following equation for gxT- 

T3 I . pole, ^2m / n^, K {v' , Q'^)aTT{v' , Q"^) , 

Re[gTT{v, Q^) - Qtt Q n = (^)^ j _ ^i ' (1) 

where is the nucleon pole (elastic) contribution, V denotes the principal value 
integral and K is the virtual photon flux factor. The lower limit of the integration 
fo is the pion-production threshold on the nucleon. A low-energy expansion gives: 

, On- 

Re[.gTT(i^, Q') - g^T%\^^ Q')] = {-^)iTT{Q')v + iTTiQ'V + 0(v^), (2) 

with a the electromagnetic fine-structure constant and M the nucleon mass. 
Itt{Q'^) is the coefficient of the 0{i^) term of the Compton amplitude. Equation 
(2) defines th e gen eralized forward spin polarizability 7tt(Q^) (or 7o(Q^) as it was 
used in Refs. ^^^201)^ Combining Eqs. ( 1) and (2), the 0{iy) term yields a sum rule 
for the generalized GDH integral ^ '^ ^ ^^ ^ 



V 



2M2 /-^"r ^ 2^ 4M2 2 / 2n1 

=-qr J [9iix,Q ) - -Q^x 92[x,Q jjdx (3) 

where x — /2Mv is the Bjorken scaling variable. As 0, the low-energy 

theorem relates 1(0) to the anomalous ma gnetic moment of the nucleon, k, and 
Eq. ([3]) becomes the original GDH sum rule^^ES 

where cr]^/2(3/2) is the total photoabsorption cross section with a projection of i 
(|) for the total spin along the direction of photon momentum, and 2aTT = 
o'i/2 — o'3/2- The 0(1/"^) term yields a sum rule for the generalized forward spin 
polarizabilitjl22I2Il 

l&aAP p 2r , ^2N 2 , ^2^lJ 

=— ge — y ^ [9i{x,Q )~-^x g2[x,Q )\^dx. (5) 

Considering the longitudinal-transverse interference amplitude gLT: with the 
same assumptions, one obtains: 

M9lt{v. Q^) - glf{v, Q^)] = (^)Q/ij.(Q2) + QS^^iQ^^ + 0{v^) (6) 

where the 0{\) term leads to a sum rule for Ilt{Q^), which relates it to the ctlt 
integral over the excitation spectrum. 
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2M2 







giix,Q^)+g2ix,Q^) 



dx. 



The 0{v^) term leads to the generahzed longitudinal-transverse polarizability 



(7) 



1 



dv 



gi{x,Q^)+g2ix,Q^) 



dx. 



(8) 



Alternatively, we can consider the covariant spin-dependent VVCS amplitudes 
5'i and 5*2, which are related to the spin- flip amplitudes grr and gLT- 



vM 



Dispersion relation with the same assumptions leads to 
where the 0(1) term leads to a sum rule for /i((5^): 



2A/2 



gi{x,Q'^)dx. 



(9) 
(10) 

(11) 
(12) 



The 0{v^) term leads to the generalized gi polarizability: 



6 



2a 



LT 



x'^gi{x, Q^)dx 

iTTiQ^) - h{Q'' 



A/2Q2V--^- ' V ^^^^ 

For 5*2, assuming a Regge behavior at — > oo given by ^2 — > i^"^ with a2 < 
^1, the unsubtracted dispersion relations for S2 and 1^52, without the elastic pole 
subtraction, lead to a "super-convergence relation" that is valid for any value of Q^, 

g2{x,Q^)dx^0, (14) 

which is the Burkhardt-Cottingham (BC) sum rultPl. It can also be written as 



2M^ 



' 52(x, Q^)dx = ^Fp{Q^) (Fd{Q'') + Fp{Q^) 



(15) 



where Fp and Fd are the Pauli and Dirac form factors for elastic e-N scattering. 
The low-energy expansion of the dispersion relation leads to 

Re [{1^82(1^, Q^)) - (i/^f°'^(i^,Q2))] 
= 2a/2(Q') - ^(/tt(Q^) - hiQ'))i^' + ^l,AQ'y + 0{v% (16) 
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where the 0{h''^) term gives the generahzed g2 polarizabiUty: 
IgAQ ) = — / ^ 92{x,Q )dx 



Jo 

= Slt{Q') - iTTiQ') + j^{Itt{Q') - h{Q')) ■ (17) 
At high the OPEj 27 | 28 | 29 l ^-j^^ VVCS amphtude leads to the twist expan- 



sion: 



Mr(Q2) , . 

2)72 (1^) 



with the coefficients fi-r related to nucleon matrix elements of operators of twist 
< r. Here twist is defined as the mass dimension minus the spin of an operator and 
/It- itself is a pertubative series in a^, the effective strong coupling constant. Note 
that the application of the OPE requires summation over all hadronic final states, 
including the elastic at a; = 1. 

The leading-twist (twist-2) component, fi2^ is determined by matrix elements 
of the axial vector operator iplfj.lbi'j summed over quark flavors, where are the 
quark field operators. It can be decomposed into flavor triplet (g^), octet (ag) and 
singlet (AS) axial charges, 

M2(Q') = (±^9A + ^as) + ^AS + OiasiQ^)), (19) 

where +(-) corresponds to proton (neutron) and the 0(qs) terms are the evo- 
lution due to the QCD radiative effects that can be calculated from perturbative 
QCD. The triplet axial charge is obtained from neutron /3-decay, while the octet 
axial charge can be extracted from hyperon weak-decay matrix elements assuming 
SU(3) flavor symmetry. Within the quark-parton model AE is the amount of nu- 
cleon spin carried by the quarks. Deep Inelastic Scattering (DIS) experiments at 
large have extracted this quantity through a global analysis of the world datclH. 

Eqs. (fTS]) and p9|l . at leading twist, lead to the Ellis- Jaffe sum rulJ^SI t^h)^ 
the assumptions of SU(3) flavor symmetry and an unpolarized strange sea. The 
difference between the proton and the neutron gives the flavor non-singlet term: 

rW) ~ ^liQ') = + o{a,) + o{i/Q% (20) 

which becomes the Bjorken sum rule at the —>■ oo limit. 

If the nucleon mass were zero, fi^ would contain only a twist-r operator. The 
non-zero nucleon mass induces contributions to jir from lower- twist operators. The 
twist-4 term contains a twist-2 contribution, 02, and a twist-3 contribution, d2, in 
addition to /2, the twist-4 component^'^ 28 29 30. 

/i4 = M2 (02 + 4rf2 + 4/2) /9. (21) 
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The twist-2 matrix element a2 is: 

a2 S^^P^P^^ = ^ E l^''iD''iD'^^,\P, S) , (22) 

where Cg is the electric charge of a quark with flavor q, S and P are the covariant 
spin and momentum vectors, D" are the covariant derivatives, and the parentheses 
{• • •} denote symmetrization of indices. The matrix element 02 is related to the 
second moment of the twist-2 part of gi: 

a2{Q^) = 2 f dxx^ gi{x, Q^) . (23) 
Jo 

Taking Eq. ([23|) as the definition of 02, it is now generalized to any Q^, including 
twist-2 and higher-twist contributions. At low Q^, the inelastic part of 02, is related 
to 7gj^ , the generalized gi polarizability: 

Note that at large Q^, the elastic contribution is negligible and Q2 becomes 02. 
The twist-3 component, 1^2, is defined by the matrix elemenlP^I^SHnEni 

d^Si^^pMp^} ^^J2{P:S\i^g gM^-^^^.lP^S) , (25) 

where g is the QCD coupling constant, F'^" = (l/2)e'^''"''Fa^, Faf3 are the gluon 
field operators, and the parentheses [• ■ •] denote antisymmetrization of indices. This 
matrix element is related to the second moments of the twist-3 part of gi and g2- 

d2(Q') =^ dx x^(2gi{x,Q^) + 3g2{x,Q^)') 

- 3j\xx'[g2{x,Q')-gY'^{x,Q')), (26) 
where gY^ is the twist-2 part of 52 as derived by Wandzura and WilczelPIl 

9r'^ix,Q')^g,ix,Q')+ fdy'-l^^ . (27) 

Jx V 

The definition of ^2 with Eq. (|26p is generalized to all Q^. At low Q^, the inelastic 
part of d2{Q'^) is related to the polarizabilities: 

At large Q^, c?2 becomes d2 since the elastic contribution becomes negligible. 
The twist-4 contribution to fi^ is defined by the matrix element 

/2 M'S'' ^ i ^ (Ar|5 Ft^-^j, ^,|7V) (29) 
9 
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where F^'^ is the dual gluon-field strength tensor. A (generalized) definition of /2 
in terms of the structure functions is 

1 /"^ 

MQ^)=2j^ dxx^{7g,ix,Q^) + 12g2ix,Q^)-9g3{x,Q^)y (30) 

where 53 is a parity-violating spin structure function which can be accessed by 

measurement of unpolarized lepton scattering off a longitudinally polarized target 

or with neutrino scattering. With only gi and 52 data available, /2 can be extracted 

through Eqs. (fT8|) and pT|) if the twist-6 or higher terms are not significant. 

The twist-3 and 4 operators describe the response of the collective color electric 

and magnetic fields to the spin of the nucleon. Expressing these matrix elements in 

terms of the components of F^'^ in the nucleon rest frame, one can relate ^2 and /2 

97 9Q '^0 

to color electric and magnetic polarizabilities. These are defined a a i^o,^u , ou i 

XE 2M^S = (iV| ja X Ea \N) , XB 2M^S = {N\ jl \N) , (31) 

where S is the nucleon spin vector, is the quark current, Ea and Ba are the 
color electric and magnetic fields, respectively. In terms of 6,2 and /2 the color 
polarizabilities can be expressed as 

XE = \ (2^2 + h), XB^\ (4^2 - h) ■ (32) 



3. Summary of previous experimental situation 

Before Jefferson Lab started running polarized beams and targets, most of the spin 
structure measurements of the nucleon we re performed at hig h-energy facilities like 
CERN (EMcjSll and SMcPEHSSIMETl)^ oggy (HERMES^^'^^ ^O") gLAC 
(E8041, E13(P, E14dSI, E14d31l, E15i'45 46^ E15547 48 ^nd E155^). The mea- 
sured gi and 32 data were suitable for an analysis in terms of perturbative QCD. 
The impetus for performing these experiments on both the proton and the neutron 
was to test the Bjorken sum rule, a fundamental sum rule of QCD. After twenty-five 
years of active investigation this goal was accomplished with a test of this sum rule 
to better than 10%. The spin structure of the nucleon was unraveled in the same 
process. Among the highlights of this effort is the determination of the total spin 
content of the nucleon due to quarks, AS (see Eq. ([T9|) ). It revealed the important 
role of the quark orbital angular momentum and gluon total angular momentum. 
The main results from this inclusive double spin asymmetry measurement program 
have led to new directions namely the quest for an experimental determination of 
the orbital angular momentum contribution^^ (^.g-, with Deep Virtual Compton 
Scattering at Jef fers on Lab and other facilities) and the gluon spin contribution 
(with COMPAS^^and RHIC-spiilSli experiments) . These efforts will be ongoing 
for the next decades. 
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4. Description of the JLab experiments 

The inclusive experiments described here took place in the fall of 1998 in JLab 
Halls aE31 and b531. The accelerator produces a CW electron beam of energy up 
to 6 GeV. In Hall A the expriment was performed using incident beam energies of 
5.06, 4.24, 3.38, 2.58 1.72 and 0.86 GeV and a beam current up to 15 /uA. In Hall B, 
beam energies of 2.56 and 4.28 GeV and currents up to 2.5 nA were used. A beam 
polarization of about 0.70, as measured by Halls A and B M0ller polarimeters, was 
obtained by illuminating a strained GaAs cathode with polarized laser light. Data 
taken with unpolarized targets show that the beam charge asymmetry and false 
asymmetries are under control. 

4.1. Hall A 

A polarized high pressure (~12 atm.) gaseous ^He targeJ^I ^^ I ^^I was used as an 
effective polarized neutron target in the experiment performed in Hall A. The target 
cell is made of two connected chambers. Rubidium atoms, confined by a thermal 
gradient to the upper spherical chamber, are polarized by optical pumping. The 
polarization is then transfered to ^He nuclei by spin-exchange collisions. The latter 
descend by diffusion and convection to the lower chamber, a 40 cm long tube of 
2 cm diameter, where they interact with the beam. The spins are held in a given 
direction by a 2.5 mT uniform field. The average target polarization, monitored 
by NMR and EPR techniques'^ ^'^ , was 0.35±0.014 and its direction could be 
oriented longitudinal or transverse to the beam direction. The measurement of cross 
sections in the two orthogonal directions allowed a direct extraction of gi^"^ and 
§2^'^, or equivalently (Jtt and ctlt, without the need of models or unpolarized data. 

The beam energy was measured with two independent devices to a few 10~^ 
precisiorP3^, while its current was monitored to better than 1%. The scattered elec- 
trons were detected by two High Resolution Spectrometers (HRS) of ~ 6 msr an- 
gular acceptance and ~ 9% relative momentum acceptance, pointing at a nominal 
15.5° angle. The high luminosity of 10^^ cm^^s"^ allowed for statistically accurate 
data at numerous HRS momentum settings that covered the [ff' ,v) plane, as is nec- 
essary to form moments at fixed . The detector package consisted of vertical drift 
chambers (for momentum analysis and vertex reconstruction) , scintillation counters 
(data acquisition trigger) and Cerenkov counters and lead-glass calorimeters (for 
particle identification (PID)). The 7r~ were sorted from e~ with an efficiency better 
than 99.9% . Both HRS spectrometers were used to double the statistics and con- 
strain the systematic uncertainties, by comparing the cross sections extracted using 
each HRS. Acceptances and optical properties were studied with a multi-foil carbon 
target and sieve-slit collimators. Carbon elastic cross sections were measured with 
2% accuracy. 

The kinematic coverage spans from the elastic reaction to slightly beyond the 
resonance region. Data from the pion threshold to W — 2 GeV were used to form 
neutron moments. Polarized '^He elastic data were used to cross-check target and 



July 29, 2008 19:27 WSPC/INSTRUCTION FILE sumrules 



Sum Rules and Moments of the Nucleon Spin Structure Functions 9 

beam polarimetry at the 4% level and absolute cross sections at the 5 — 7% level. 
Electromagnetic radiative corrections were performed using the method of Mo and 
Tsai^ for external radiative corrections, while internal radiative corrections were 
evaluated with a modified version of the POLRAD code^ including new quasi- 
elastic and resonance data. The relative uncertainty on the radiative corrections 
was typically less than 20%. Other systematic uncertainties came from absolute 
cross-section measurements (5%), beam and target polarizations (both relative 4%). 

The spin structure functions g" and g2 are extracted using polarized cross- 
section differences in which contributions from unpolarized materials such as target 
windows, nitrogen or the two protons of the "^He nucleus cancel. Corrections for the 
two protons in "^He are still needed since they are slightly polarized due to the D state 
(~ 8%) and S' state (~ 1.5%) of the ^He wave function^. Corrections for binding 
and Fermi motion are also applied via a PWIA-based model^. The uncertainty on 
this correction varies from 10% at the lowest to < 5% at larger Q^. To form 
the neutron moments, the integrands (e.g. axT or gi) are needed at constant Q^. 
Six equally spaced values of were chosen in the range 0.1 < < 0.9 GeV^ 
and the integrands were determined from the measured points by interpolation. 
To complete the moments for the unmeasured high-energy region, the Bianchi and 
Thomas parameterizatiorP^ was used for 4 < < 1000 GeV^ and a Regge-type 
parameterization was used for > 1000 GeV^. 



4.2. Hall B 

Polarized solid ^^NHs and ^^NDa target using dynamic nuclear polarizatiorl^lEl] 
were used in Hall B. Polarizations varied from 70% to 40% for NH3 and 25% to 
10% for ND3. The target material was cooled to 1.2 K and the spins were oriented 
by a longitudinal 5 T field. 

The beam was rastered to evenly spread depolarization effects. The beam energy 
was obtained from Hall A measurements. Its current was monitored by a Faraday 
cup downstream of the target. The luminosity reached 10'^^cm~^s~^. The CEBAF 
Large Acceptance Spectrometer (CLAS) has a large angular (2.57r sr) and momen- 
tum acceptance, which is ideal for studying the evolution of moments. CLAS 
contains six superconducting coils that produce a toroidal magnetic field. Particles 
are characterized by three layers of drift chambers (tracking, momentum and charge 
information), a layer of scintillator paddles (time of flight), and Cerenkov counters 
supplemented by electromagnetic calorimeters (PID). After PID, the remaining tt~ 
contamination was less than 1% of the electron rate. To cover lower angles and 
momenta, the target was shifted 55 cm upstream of its nominal location. The polar 
angle ranged from 8° to 50°. The azimuthal angular acceptance was about 80%. 
The lowest momentum accepted was 0.35 GeV for the proton run. The deuteron 
data were taken with in-bending torus field only, which yields a minimum angle of 
14° and a minimum momentum of 0.5 GeV. 

Since cross sections are more difficult to measure in a large-acceptance spec- 
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trometer and F2{x,Q^) and R{x,Q^) are relativel y w ell it is easier to 

extract gi{x,Q'^) using asymmetry measurements^. The physics asymmetry A\\ 
is related to the measured one: Araw — A\\{Pi)PtDF) / Cn ■ CAr=0.98 is the effect 
of the polarized valence proton of ^^N , Pb{Pt) is the beam (target) polarization. 
DF 0.15(0.2) for NHg (ND3) is the dilution factor from ^^N, ''He and windows 
on the beam path. DF was estimated using data from a carbon target. Parameteri- 
zations of proton and neutron cross sections were used to account for the difference 
between ^^N and ^^C nuclei. Acceptance and detector inefficiency cancels out to 
first order in the asymmetry ratio. Acceptance cuts were chosen such that second- 
order corrections remain negligible. The radiative corrections to yl|| were done using 
RCSLACP0#6 based on the same formalism as POLRAD. The cross sections are 
parameterized using the NMC fit of the unpolarized structure-function datcP3 and 
the polarized data from SLAG, CERN and HERME^. The resonance contribution 
is estimated by the MAIeJI^ and the "AO"^^ models. The radiative corrections 
account for 20% of the total systematic uncertainty on Fi at low and 5% at 
larger Q^. 

The same model was used to extract Ai and gi from A\\. Q^) RP{x, Q^) 

were measured at SLACP and JLall^. The DIS part of 52 is calculated using 
Eq. ([27|) while its resonance part is obtained from the MAID and "AO" models. 
The model-dependence dominates the systematic uncertainty on Fi, contributing 
to about 75% at low and 50% at larger for the proton and 50% in average 
for the deuteron. The model was also used to estimate the unmeasured high-energy 
part of F^''^. Nuclear corrections were applied to Tf to extract F" by accounting 
for the deuteron D-state: F" = 2F^/(1 — 1.5a;D) — Ff , with the D-state probabihty 
UJD ^ 0.05. 

5. Recent results from Jefferson Lab 

5.1. Results of the generalized GDH sum for the neutron 

Fig. [1] shows the extended GDH integrals I{Q^) (open circles) for the neutron, 
which were extracted from JLab experiment E94-010 (Hall A)P, from pion threshold 
to W = 2.0 GeV. The uncertainties, when visible, represent statistics only; the 
systematics are shown by the grey band. The solid squares include an estimate of 
the unmeasured high-energy part. The corresponding uncertainty is included in the 
systematic uncertainty band. These data indicate a smooth variation of I{Q^) to 
increasingly negative values as varies from 0.9 GeV^ towards zero. The data 
are more negative than the MAID model calculatiorP^. The calculation includes 
contributions to I{Q^) for W < 2 GeV, and should thus be compared with the open 
circles. At high Q^, the data approach the HERMES^^ neutron results (extracted 
from '^He), which spans the range 1.28 GcV^ <Q^ < 7.25 GeV^ but includes only 
the DIS part of the GDH integral. The GDH sum rule prediction, /(O) = -232.8 ^b, 
is indicated on Fig. [1] along with extensions to > using two xPT calculations, 
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Fig. 1. Results for I(Q^) are compared with xPT calc ulat ions of ref.EU (dotted line) and ref. 1221 
(dot-dash ed line). The MAID model calculation of ref . 113, is shown with a solid line. Data from 
HERMESi^SI are also shown and a log scale is used for > 1 GeV^. 



one using the Heavy Baryon approximation (HB^PT) (dotted line) and the other 
Relativistic Baryon xPT (RB^PT^^ (dot-dashed line). Shown with a grey band is 
RBxPT including resonance efFectJ^, which have an associated large uncertainty 
due to the resonance parameters used. The uncertain ty c an be seen to encompass 
both the lowest point and the calculation of Ref. ED. At = 0.3 GeV^, the 
prediction of Ref. is much more negative than the data. Further calculations as 
well as further measurementtQ will help clarify the situation. 



5.2. First moments of gi, g2 and the Bjorken sum 

Results on f i((5^) at low to moderate are available on the protorP and the 
neutron (using '^H^ and deuterorP). The integrals f^, F" , F" and F?"" are pre- 
sented in Fig. [2] along with the data (open symbols) from and 
HERMEdSSlM The integral (f J ) was formed by summing gl {g"^) from pion 
threshold up to VF = 2.6 (2.0) GeV. Nuclear corrections were applied to the deuteron 
and ^He data to extract F" as described in section SI These results are shown by 
the open symbols. The sums including an estimation of the DIS contributions are 
shown by the solid symbols. The inner uncertainty indicates the statistical uncer- 
tainty while the outer one is the quadratic sum of the statistical and systematic 
uncertainties. 

Data from the proton and neutron can be combined to form the flavor non- 
singlet moment T\{Q^) — F"((5^) predicted in the Bjorken limit by the Bjorken sum 
rulP. It is generalized at large to moderate Q using OPE that yields a relatively 
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Fig. 2. Results for (top left), T\ (top right) and the Bjorken sum (bottom left). The open 
circles (diamonds) represent the SLAC (HERMES) data. The open squares (stars) are the JLab 
CLAS (Hall A) data integrated over the resonance region, while the corresponding full symbols 
represent the full integrals. The slopes at Q^=0 predicted by the GDH sum rule arc given by the 
dotted lines. The MAID model predictions that includes only resonance contributions are shown 
by the plain lines while the dashed (dot-dashed) lines are the predictions from the Soffer-Teryaev 
(Burkert-Ioffe) model. The leading twist (LT) Q^-evolution of the moments is given by the grey 
band. The plain lines (bands) at low are the xPT predictions by Ji et al. (Bernard et al.). 
is shown on the bottom right panel. 

simple leading-twist evolution in which only non-singlet coefficients survive: 
f « - = ||1 - - 3.58 (^)^ - 20.21 (!^)' + ...1 + 0(±). (33) 

The Bjorken sum rule was verified to a level better than 10% at = 5 GeV^. 

Considering a flavor non-singlet object results in additional simplifications: the 
A resonance contri but ion mostly cancels out so that xPT predictions are expected 
to be more reliabltjSI], with possibly an extended domain of validity in Q^. Fur- 
thermore, a flavor non-singlet moment is a convenient quantity for Lattice QCD 
calculations since disconnected diagrams, which are difficult to calculate, do not 
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contribute. Hence, the Bjorken sum is particularly suited to provide benchmark 
tests for the three theoretical frameworks used to study the transition from hadronic 
to partonic degrees of freedom. 

To form r^~", the proton and neutron data were re-analyzed at the same 
bins. The unmeasured low- a; part of the integral was re-evaluated for p, d and '^He 
data sets using a prescription described in section [4. II with an additional constraint 
imposed by the Bjorken sum rule at = 5 GeV^. 

At Q^=0, the GDH sum rule predicts the slopes of moments (dotted lines). 
The deviation from the slopes at low can be calculated with xPT. We show 
calculations by Ji et o/!^ using the HBxPT and by Bernard et al. witli^ and 
withoulP2l including vector mesons and A degrees of freedom (the band shows a 
range of results due to the uncertainty in the 7NA form factor.) We do not show 
this result for the Bjorken sum because of the unknown correlation in the systematic 
uncertainties of the individual Ti and Ti . For the neutron this calculation overlaps 
with a data point at — 0.1 GcV^ while none of the xPT calculations agree with 
the proton data, of which the lowest point is 0.17 GeV^. The calculation of Ji 
et al. agrees with the lowest data for the Bjorken sum. At moderate and large 

we show the MAID calculatiorff^, to be compared with the data summed up to 
W — 2 GeV. The model disagrees for Ti, and Ti , with the wrong sign of the 
slope at the photon point for " due to its underestimate of the GDH integral for 
the neutron at = 0. The disagreement indicates that other contributions, not 
included in the model, are important. 

The other calculations shown are by Soffer and Teraye and by Burkert and 
loff^^. Those should be compared to the full experimental sum. The Soffer and 
Terayev model relies on the smooth Q^-behavior of the It — Ii + I2 integral to 
interpolate it between its known — and large value. Fi is extracted using 
the BC sum rule prediction on F2, Eq. (fT5|l . Power corrections and pQCD radiative 
corrections were recently added to the model to resolve the initial discrepancy with 
the JLab data. The Burkert and loffe model relies on a parameterization of reso- 
nances and non-resonance background^ supplemented by a DIS parameterization 
based on vector meson dominanc4^. Both models agree well with the data. 

The leading-twist pQCD evolution is shown by the grey band. For F^*' ' we used 
AE=0.15(Ai;=0.35), a8=0.579 and 5a=l-267 to anchor the evolution (see Eqs. (HH) 
and (|33p ). The Bjorken sum rule sets the absolute scale for F^ ".In all cases, the 
leading twist tracks the data down to surprisingly low Q^, which indicates an overall 
suppression of higher-twist effects. This is quantatively discussed in section [5.41 

The capability of transverse polarization of the Hall A '^He target allows precise 
measuremetns of 32- The integral F2 is extracted and plotted in the bottom-right 
panel of Fig. [2] in the measured region (open stars) and after adding the elastic con- 
tribution (solid stars )P^. The solid triangles correspond to the results obtained after 
adding to the solid stars an estimated DIS contribution assuming 52 — ■ The 
MAID estimate (solid line) agrees well with the resonance data. The positive light 
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grey band corresponds to the total experimental systematic errors. The negative 
dark band is the estimate of the systematic error for the low- a; extrapolation. The 
results (solid stars) are consistent with the BC sum rule to within 2 standard de- 
viations over the measured range. The SLAC E155x collaboratiorPS previously 
reported a neutron result at high (solid circle), where the elastic contribution is 
negligible; the result is consistent with zero but with a rather large error bar. On 
the other hand, the proton result was reported to deviate from the BC sum rule by 
3 standard deviations. The quoted uncertainty in this case is 3 times smaller than 
that for the neutron, but still large. 

5.3. Spin Polarizabilities: 70, S^t o-nd d-2, for the neutron 

The generalized spin polarizabilities provide benchmark tests of xPT calculations 
at low and of Lattice QCD calculations at moderate to high Q^. Since the gener- 
alized polarizabilities have an extra weighting compared to the first moments 
(GDH sum or Ilt)i these integrals have less contribution from the large-i' region 
and converge much faster, which minimizes the uncertainty due to the unmeasured 
region at large v. 

At low Q^, t he generalized polarizabilities have been evaluated with xPT 
calculation9^^21221. One issue in the xPT calculations is how to properly include the 
nucleon resonance contributions, especia lly the A resonance, which usually domi- 
nates. As was pointed out in Ref. I^^ l ^^i while 70 is sensitive to resonances, 5lt is 
insensitive to the A resonance. Measurements of the generalized spin polarizabil- 
ities are an important step in understanding the dynamics of QCD in the chiral 
perturbation region. 

The first results for the neutron generalized forward spin polarizabilities 7o(Q^) 
and 5lt{Q^) were obtained at Jefferson Lab Hall A'' over a range from 0.1 to 
0.9 GeV^. We will first focus on the low region where the comparison with xPT 
calculations is meaningful, and then show the complete data set. The results for 
7o(Q^) are shown in the top-left panel of Fig. [3]for the two lowest values of 0.10 
and 0.26 GeV^. The statistical uncertainties are smaller than the size of the symbols. 
The data are compared with a next-to-leading order, O(p^), HB^PT calculatiorP^l, 
a next-to-leading-order RB^PT calculation' , and the same calculation explicitly 

no 

including both the A resonance and vector meson contributions . Predictions from 
the MAID modeff^ 

are also shown. At the lowest point, the RB^PT calculation 
including the resonance contributions is in good agreement with the experimental 
result. For the HBxPT calculation without explicit resonance contributions, dis- 
crepancies are large even at — 0.1 GeV^. This might indicate the significance of 
the resonance contributions or a problem with the heavy baryon approximation at 
this . The higher data point is in good agreement with the MAID prediction, 
but the lowest data point at = 0.1 GeV^ is significantly lower, consistent with 
what was observed for the generalized GDH integral result (Section [HH]). Since 5lt 
is insensitive to the dominating A resonance contribution, it was believed that 5lt 



July 29, 2008 19:27 WSPC/INSTRUCTION FILE sumrules 



Sum Rules and Moments of the Nucleon Spin Structure Functions 15 



0.0 



o 



-5.0 
3.0 K 
= 2.0 




'H'MOIO 

MAIDe: 

Kao el ai 0(p VO(//) 

■ Bernard el al. 

■■ Bernard e/(;/.(VM + A) / 




-0.01 



0.002 



1 






1 ■ E94010 




3 □ SLAC 1 




8 O Lattice QCD 




i MAID 




J Kao el al. 0(p)+0(p') 




J Bemai-d et al. 

f ^ Bemai-d ffa/. (VM+A) 










1. . .. 



O.I 0.2 

e=(GeVT 



10 



G" (GeV") 



Fig. 3. Results for the neutron spin polarizabilities 70 (top panels) and Si^t (bottom panels). 
Solid squares are the results with statistical uncertainties. The light bands are the systematic 
uncertainties. The dashed curves are the HBxPT calculatiorT^ 
dark bands are the RBxPT calculatioriH], without and witl 



The dot-dashed curves and the 
ZSI the A and vector meson contri- 



butions, 



respectively. Solid curves are from the MAID modeffZl, The right panels are with 



weighting. The open squares are the SLAC data^S and the open diamonds are the Lattice QCD 
calculation; 



80 



should be nior e suita ble than 70 to serve as a testing ground for the chiral dy- 
namics of QCEl22Eni. The bottom-left panel of Fig. [3] shows Slt compared to xPT 
calculations and the MAID predictions. It is surprising to see that the data are in 
significant disagreement with the xPT calculations even at the lowest Q^, 0.1 GeV^. 
This presents a significant challenge to the present theoretical understanding. The 
MAID predictions are in good agreement with the results. 

The right panels of Fig. [3] show 70 and 6lt multiplied by along with the 
MAID and xPT calculations. Also shown are the world data^^ and a quenched 
Lattice QCD calculatiorPSl^ both at = 5 GcV^. It is expected that a t large Q^, 
the Q^-weighted spin polarizabilities become independent of (scaling|2ni21]^ ^^^^ 
the DIS Wandzura-Wilczek relational! leads to a relation between 70 and S^t'- 



Slt{Q^ 



1 



70 (Q' 



(34) 



The results show that scaling is not observed for < 1 GeV^. For the higher 
moments scaling is expected to start at a higher than for the first moments, 
for which scaling was observed to start around of 1 GeV^, where higher- twist 
effects become insignificant. Again, both results are in good agreement with the 
MAID model. 

Another combination of the second moments, d2{Q^), provides an efficient way 
to study the high behavior of the nucleon spin structure, since it is a matrix 
element, related to the color polarizabilities and can be calculated from Lattice 
QCD. It also provides a means to study the transition from high to low Q^. In 
Fig. m d2{Q^) is shown. The experimental results are the solid circles. The grey 
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band represents the systematic uncertainty. The world neutron results froni^^open 
square) and from JLab E99-ll'i^S2! (solid square) are also shown. The solid line is 
the MAID calculation containing only the resonance contribution. At low the 
HBxPT calculatioiPS' (dashed line) is shown. The RBxP"^ with or without the 
vector mesons and the A contributionJ^ are very close to the HBxPT curve at this 
scale, and are not shown on the figure for clarity. The Lattice QCD predictiorPSl at 
q2 _ 5 GeV^ is negative but close to zero. We note that all models (not shown at 
this scale) predict a negative or zero value at large Q^. At moderate Q^, our data 
show that d2 is positive and decreases with . 
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Fig. 4. The left panel shows the JLab results of (i2 for the neutron along with the world data at 
high Q^, Lattice QCD, MAID model and H ByP T calculations. The right plot compares w ith 
calculations of Stein et a/.i2Sl, MIT bag modeSHI, QCD sum ruleJSSI and an instanton modelSSI. 



5.4. Higher Twist extractions and 

The higher-twist contribution to Fi can be obtained by a fit with an OPE series, 
Eq. (|18p . truncated to an order appropriate for the precision of the data. The goal 
is to determine the twist-4 matrix element /2. Once /i4 is obtained, /2 is extracted 
by subtracting the leading- twist contributions of 02 and d2 following Eq. (HI]). To 
have an idea how the higher-twist terms (twist-8 and above) affect the twist-4 term 
extraction, it is necessary to study the convergence of the expansion and to choose 
the range in a region where /is term is not significant. This study is made possible 
only with the availability of the new low-Q^ data from JLab. 

Higher-twist an alyses have been performed on the proton , the neutron and 
the Bjorken su nP. An earlier proton analysis is aval lableP but will not be presented 
here since it uses a different procedure. Fi at moderate is obtained as described 
in section H) For consistency, the unmeasured low- a; parts of the JLab F^ and of 
the world data on Fi were re-evaluated using the same prescription previously 
used for F" and F^^". The elastic contribution, negligible above of 2 GeV^ but 
significant (especially for the proton) at lower values of Q^, was added using the 
parametrization of Ref . The leading-twist term ^2 was determined by fitting the 
data at > 5 GeV^ assuming that higher twists in this region are negligible. 
A value of qa = 1-270 ± 0.045 was obtained for the Bjorken sum. Using the proton 
(neutron) data alone, and with input of qa from the neutron beta decay and ag from 
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the hyperon decay (assumed SU(3) flavor symmetry), we obtained AE = 0.15±0.07 
for the proton and AS = 0.35 ± 0.08 for the neutron. We note that there is a 
significant difference between AS determined from the proton and from the neutron 
data. This is the main reason why the extracted 114 and /2 from the Bjorken sum is 
different compared to the difference of those extracted individually from the proton 
and neutron, since the Bjorken sum does not need the assumption of SU(3) flavor 
symmetry and AS was cancelled. 

The fit results using an expansion up to (l/Q^) in determining ^4 are summa- 
rized in Table[TJ In order to extract /2, shown in Table[2l the target-mass corrections 
a2 were evaluated using the Blumlein-Boettcher world data parametrizatiorP^l for 
the proton and a world fit to the ne utron data, which includes the recent high pre- 
cision neutron results at large af^. The ^2 values used are from SLAC EISSjP^ 
(proton) and JLab E99-117 (neutron). 



Table 1. Results of fia and fig at = 1 GeV^ for proton, neutron and p — n. The uncertainties 
are first statistical then systematic. 



Target 


(GeV^) 








proton 
neutron 

p — n 


0.6-11.0 
0.5-11.0 
0.5-11.0 


-0.065±0.012 ± 0.048 
0.019±0.002 ± 0.024 
-0.060±0.045 ± 0.018 


0.143±0. 021 ±0.056 
-0.019±0.002 ±0.017 
0.086±0.077 ± 0.032 


-0.026±0.008± 0.016 

0.00±0.00 ± 0.03 
0.011±0.031 ±0.019 



Table 2. Results of /2, XE and XB at = 1 GeV'^ for proton, neutron and p — n. The 
uncertainties are first statistical then systematic. 



Target 


/2 


XE 


XB 


P 
n 

p — n 


-0.160 ±0.028 ±0.109 
0.034 ±0.005 ± 0.043 
-0.136 ±0.102 ±0.039 


-0.082 ±0.016 ±0.071 
0.031 ±0.005± 0.028 
-0.100 ±0.068 ±0.028 


0.056 ±0.008 ± 0.036 
-0.003 ± 0.004 ± 0.014 
0.036 ±0.034 ±0.017 



The Fi were fit, varying the minimum threshold to study the convergence 
of the OPE series. The extracted quantities have large uncertainties (dominated by 
the systematic uncertainty) but are stable with respect to the minimal threshold 
when it was below 1 GeV^. The results do not vary significantly when the /ig term 
is added, which justifies a posteriori the use of the truncated OPE series in the 
chosen range. In the proton case, the elastic contribution makes a significant 
contribution to the fie term at low but this does not invalidate a priori the 
validity of the series since the elastic contributes mainly to /ig and /is remains small 
compared to fi4. We notice the alternation of signs between the coefficients. This 
leads to a partial suppression of the high er-twist effects and may be a reason for 
quark-hadron duality in the spin sector^. We also note that the sign alternation is 
opposite for the proton and neutron. The results are compared to theoretical calcu- 
lations in Fig.m Following Eq. ([5^ . the electric and magnetic color polarizabilities 
were determined. Overall, the values, given in Tabled are small, and we observe a 
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sign change in the electric color polarizability between the proton and the neutron. 
We also expect a sign change in the color magnetic polarizability. However, with 
the large uncertainty and the small negative value of the neutron xb, it is difficult 
to confirm this expectation. 

6. Conclusion 

A large body of nucleon spin-dependent cross-section and asymmetry data have 
been collected at low to moderate in the resonance region. These data have 
been used to evaluate the evolution of moments of the nucleon spin structure 
functions gi and (/2, including the GDH integral, the Bjorken sum, the BC sum 
and the spin polarizabilities, and to extract higher-twist contributions. The latter 
provided access to the color polarizabilities in the nucleon. 

At close to zero, available next-to-leading order xPT calculations were tested 
against the data and found to be in reasonable agreement for — 0.1 GeV^ 
for the GDH integral I{Q^), Ti{Q'^) and the forward spin polarizability "foiQ^)- 
Above =0.1 GcV^ a significant difference between the calculation and the data 
is observed, pointing to the limit of applicability of xPT as becomes larger. 
Although it was expected that the xPT calculation of 6lt would off'er a faster 
convergence because of the absence of the A contribution, the experimental data 
show otherwise. None of the available xPT calculations can reproduce Slt at 
= 0.1 GeV^. This discrepancy presents a significant challenge to our theoretical 
understanding at its present level of approximations and might indicate that higher 
order calculations arc necessary. 

Overall, the trend of the data is well described by phonemenological models. 
The dramatic evolution of Ig dh from high to low was observed as predicted 
by these models for both the proton and the neutron. This behavior is mainly 
determined by the relative strength and sign of the A resonance compared to that 
of higher energy resonances and deep inelastic processes. This also shows that the 
current level of phenomenological understanding of the resonance spin structure 
using these moments as observables is reasonable. 

The neutron BC sum rule is observed to be verified within uncertainties in the 
low range due to a cancellation between the resonance and the elastic contri- 
butions. The BC sum rule is expected to be valid at all Q^. This test validates the 
assumptions going ino the BC sum rule, which provides confidence in sum rules 
with similar assumptions. 

In the region above 0.5 GeV^ the first moment of gi for the proton, neutron 
and the proton-neutron difference were re-evaluated using the world data and the 
same extrapolation method of the unmeasured regions for consistency. Then, in the 
framework of the OPE, the higher twist /2 and ^2 matrix elements were extracted. 
The low data allowed us to gauge the convergence of the expansion used in this 
analysis. The extracted higher- twist (twist-4 and above) effects arc not significant 
for above 1 GeV^. This fact may be related to the observation that quark-hadron 
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duality works reasonably well at above 1 GeV^. 

Finally, the proton and neutron electric and magnetic color polarizabilities were 
determined by combining the twist-4 matrix element /2 and the twist-3 matrix 
element ^2 from the world data. Our findings show a small and slightly positive 
value of xe and a value of xb close to zero for the neutron, while the proton has 
slightly larger values for both xe and xb but with opposite signs. 

Overall, the new JLab data have provided valuable information on the transition 
between the non-perturbative to the perturbative regime of QCD. They form a 
precise data set for twist expansio n an alysis and a check of xPT calculations. 

Future precision measurement^!^ of the gi and g2 structure functions at w 
1 — 4 GeV^ will reduce the uncertainty in the extracted higher-twist coefficients and 
provide a benchmark test of Lattice QCD. 
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